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We investigate cyclic and singularity-free evolutions in a universe governed by Lagrange-multiplier 
modified gravity, either in scalar-field cosmology, as well as in f[R) one. In the scalar case, cyclicity 
can be induced by a suitably reconstructed simple potential, and the matter content of the universe 
can be successfully incorporated. In the case of /(_R)-gravity, cyclicity can be induced by a suitable 
reconstructed second function /2 (7?) of a very simple form, however the matter evolution cannot be 
analytically handled. Furthermore, we study the evolution of cosmological perturbations for the two 
scenarios. For the scalar case the system possesses no wavelike modes due to a dust-like sound speed, 
while for the f[R) case there exist an oscillation mode of perturbations which indicates a dynamical 
degree of freedom. Both scenarios allow for stable parameter spaces of cosmological perturbations 
through the bouncing point. 
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I. INTRODUCTION 

Inflation is now considered to be a crucial part of the 
cosmological history of the universe [l|, however the so 
called "standard model of the universe" still faces the 
problem of the initial singularity. Such a singularity 
is unavoidable if inflation is realized using a scalar field 
while the background spacetime is described by the stan- 
dard Einstein action Q. As a consequence, there has 
been a lot of effort in resolving this problem through 
quantum gravity effects or effective field theory tech- 
niques. 

A potential solution to the cosmological singularity 
problem may be provided by non-singular bouncing cos- 
mologies [3|- Such scenarios have been constructed 
through various approaches to modified gravity, such as 
the Pre-Big-Bang jj| and the Ekpyrotic [5|_models, grav- 
ity actions with higher order corrections [6|, braneworld 
scenarios [m, non-relativistic gravity [8|, loop quantum 
cosmology [9| or in the frame of a closed universe |10| . 
Non-singular bounces may be alternatively investigated 
using effective field theory techniques, introdu cing mat- 
ter fields violating the null ener gy c ondition (ll| lead- 
ing to observable predictions too |12| . or introduce non- 
conventional mixing terms |13| . The extension of all the 
above bouncing scenarios is the (old) paradigm of cyclic 
cosmology |14| , in which the universe experiences the pe- 
riodic sequence of contractions and expansions, which has 
been rewaked the last years [1^ since it brings different 
insights for the origin of the observable universe |16l4l8| 
(see [l^ for a review). 

On the other hand, the interesting idea of modifying 
gravity using Lagrange multipliers, first introduced in 



[3| with bouncing solution studied in [20|, was recently 
gained new interest, since its extended version may re- 
duce irrelevant degrees of freedom in modified gravity 
models |2ll - (23| . In particular, using two scalar fields, one 
of which is the Lagrange multiplier, leads to a constraint 
of special form on the second scalar field, and as a result 
the whole system contains a single dynamical degree of 
freedom. 

In the present work we are interested in constructing 
a scenario of cyclic cosmology in universe governed by 
modified gravity with Lagrange multipliers. We do that 
both in the case of scalar cosmology, as well as in the 
case of /(i?)-gravity, both in their Lagrange-multiplier 
modified versions. As we show, the realization of cyclicity 
and the avoidance of singularities is straightforward. 

This paper is organized as follows. In section |TT] we 
present modified gravity with Lagrange multipliers, sep- 
arately for scalar cosmology (subsection III A|l and for 
/(i?)-gravity (subsection IIIB[) . In section IIIII we con- 
struct the scenario of cyclicity realization, reconstructing 
the required potential (subsection IIII A| ) and the form of 
f{R) modification (subsection llll Bp . In section lTVl we in- 
vestigate the stability of the cosmological perturbations 
of the two scenarios. Finally, section |V] is devoted to the 
summary of our results. 



II. MODIFIED GRAVITY WITH LAGRANGE 
MULTIPLIERS 



Let us present the cosmological scenarios with La- 
grange multipliers. In order to be general and complete, 
we present both the scalar case, as well as the f{R)- 
gravity one. Throughout the work we consider a fiat 
Friedmann-Robertson- Walker geometry with metric 
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ds^ 



-dt^ + a^{t)dyp' 



(1) 



where a is the scale factor, although we could straight- 
forwardly generalize our results to the non-flat case too. 



A. Scalar cosmology with Lagrange multiplier 

In this subsection, we consider a scenario with two 
scalars, namely and A, where the second scalar is a 
Lagrange multiplier which constrains the field equation 
of the first one. The action reads |2ll u2l : 



S = 



a xJ—q { — d. 



(i)^^'4> - V{(t>) 



-A 



-d^cj^d^^c^ 



u{4>) 



+ C,n \ , (2) 



where A is the Lagrange multiplier field. Furthermore, 
V{4)) and U{(t>) are potentials of 0, while the function 
a;(0) and especially its sign determines the nature of the 
scalar field (/), that is if it is canonical or phantom. Fi- 
nally, the term £,„ accounts for the matter content of the 
universe. 

Variation over the A-field leads to 



o = ^-u{4>). 



(3) 



As we observe, this equation acts as a constraint for the 
other scalar field </>, and it is the cause of the significant 
cosmological implications of such a construction. Now, 
the Friedmann equations straightforwardly write 



:H^ = 



uj{(t>) + \ 



02 + y(0) + \u{(i>) 



[w(0) + 2A] U{(j)) + 1/(0) 



(4) 



2H + iH^] = 



;(</.) + A' 



02 - y(0) - A[/(0) 



Lu{<p)Ui4>)^V{^)+p„ 



(5) 



where we have defined the Hubble parameter as iJ = -, 
and we have also made use of the constraint ([3]). Ad- 
ditionally, in these expressions pm and pm stand respec- 
tively for the energy density and pressure of matter, with 
equation-of-state parameter Wm = Pm/Pm- Observing 
the above equations, we can straightforwardly define the 
energy density and pressure for the dark-energy sector as 

Pde=[c^(0)+2A][/(0)+F(0), (6) 

Pd, = c^(0)[/(0) - F(0), (7) 

and thus the dark energy equation-of-state parameter will 
be 



Wde 



C^(0)t/(0) - F(0) 

[a;(0) + 2A]C/(0) + F(0)' 



(8) 



Let us now explore some general features of the sce- 
nario at hand following [22|. First of all the constraint 



equation ^ can be integrated for positive U{(t)), leading 
to 



t = ± 



</- 



./WW) 



(9) 



Inverting this relation with respect to 0, one can find 
the explicit f-dependence of 0, that is the corresponding 
0(i). Thus, substituting the expression of 0(f) into ([5]), 
we obtain a differential equation for H[t): 

2H{t) + mit)? - -«' {^m))um)) 

-1/(0(t))+w™p„(t)}, (10) 

the solution of which determines completely the cosmo- 
logical evolution. Finally, substituting 0(i) and H{t) into 
(JD) , we can find the t-dependence of the Lagrange multi- 
plier field A: 



m = 



2t/(0(t)) \k 



^H{tf-V{4>{i))-p„,{t) 



w(0(t)) 



(11) 



However, one could alternatively follow the inverse pro- 
cedure, that is to first determine the specific behavior of 
Hit) and try to reconstruct the corresponding potential 
y(0), which is responsible for the 0-evolution that leads 
to such a H{t). In particular, he first chooses a suit- 
able t/(0) which will lead to an easy integration of ©, 
allowing for an interchangeable use of and t through 
0(f) <-!• t(0). Therefore, the second Friedmann equation 
([5]) gives straightforwardly 



F(0) 



1 



2iJ(i(0)) + 3iJ(<(0))2 

+w(0)f/(0)+u;,„p™(f(0)), (12) 



where a;(0) can still be arbitrary. In summary, such a 
potential ^(0) produces a cosmological evolution with 
the chosen H{t). 

Finally, let us make a comment here concerning the na- 
ture of models with a second, Lagrange-multiplier field, 
since it is obvious that the constraint ^ changes com- 
pletely the dynamics, comparing to the conventional 
models. In principle, one expects to have a propagation 
mode for each new field. However, due to the constraint 
([5]) and the form in which the A-field appears in the ac- 
tion, the propagating modes of and A do not appear, 
and the system is driven by a system of two first-order 
ordinary differential equations, one for each field. As a 
consequence, there are no propagating wave-like degrees 
of freedom and the sound speed for perturbations is ex- 
actly zero irrespective of the background solution J2l|. 



B. /(i?)-gravity with Lagrange multiplier 

In this subsection we present /(i?)-gravity with La- 
grange multipliers, following [22|. In this case, we start 



by a conventional /(i?)- gravity, and we add a scalar field 
A which is a Lagrange multiplier leading to a constraint. 
In particular, the action reads: 



S 



d^Xy 



-.9 ./i(i?)-A 






/2(i?) 



(13) 
where fi{R) and f2{R) are two independent functions of 
the Ricci scalar R. Note that in the above action we have 
not included the matter content of the universe, since this 
would significantly modify the multiplying terms of A, 
making the subsequent reconstruction procedures tech- 



nically very complicated. These subtleties are caused by 
the fact that the Lagrange multiplier field propagates in 
this case, which is a disadvantage of the present scenario, 
contrary to the scalar cosmology of the previous subsec- 
tion. 

Variation over A leads to the constraint 



'2^ 



f2{R). 



(14) 



Additionally, varying over the metric g^^i, and keeping the 
(0, 0)-component we obtain 



J 



0=-2/i(i?) 



18A(i? + AHHY + \'i[H + H 



-^H-\\h^^-\h^^ 



S-- 



dR 



(15) 



where ^r denotes the derivative with respect to the Ricci 
scalar. 

For J^iR) > 0, the constraint ([T^ can be solved as 



t 



dR 



^/vW 



(16) 



and inverting this relation with respect to R one can 
obtain the explicit i-dependence of i?, that is the corre- 
sponding R{i). On the other hand, the Ricci scalar is 
given from R = QH + 12 H^ . Thus, inserting R{t) in this 
relation one obtains a differential equation in terms of 
H{t), namely 



%H{t) + l2[H{t)]'^ == R{t), 



(17) 



the solution of which determines completely the cosmo- 
logical behavior. Finally, using the obtained H{t) and 
i?(t), equation ([T5]) becomes a differential equation for 
the Lagrange multiplier field A(t) which can be solved. 

However, one could alternatively follow the inverse pro- 
cedure, that is to first determine the specific behavior 
of H{t) and try to reconstruct the corresponding f2{R), 
which is responsible for the _R-evolution that leads to such 
an H{t). In particular, with a known H{t) ([TT]) provides 
immediately R{t), which can be inverted, giving t = t{R). 
Thus, using the constraint ^^, the explicit form of fiiR) 
is found to be 



«-'^K^ 



(18) 



t=t(R) 



It is interesting to mention that in the above discussion 
fi{R) is arbitrary. That is, while in conventional f{R)- 
gravity the cosmological behavior is determined com- 
pletely by fi{R), in Lagrange-multiplier modified f{R)- 
gravity the dynamics is determined completely by f2{R). 
In this case fi{R) would become relevant only in the 
presence of matter, and its effect on Newton's law |23 |. 



III. COSMOLOGICAL BOUNCE AND CYCLIC 
COSMOLOGY 



Having presented the cosmological models with La- 
grange multipliers, both in the scalar, as well as in the 
/(i?)-gravity case, in this section we are interested in in- 
vestigating the bounce and cyclic solutions. 

In principle, in order to acquire a cosmological bounce, 
one has to have a contracting phase {H < 0), followed 
by an expanding phase {H > 0), while at the bounce 
point we have H = 0. In this whole procedure the time- 
derivative of the Hubble parameter must be positive, that 
is iJ > 0. On the other hand, in order for a cosmological 
turnaround to be realized, one has to have an expanding 
phase [H > 0) followed by a contracting phase {H < 0), 
while at the turnaround point we have H — 0, and in 
this whole procedure the time-derivative of the Hubble 
parameter must be negative, that is if < 0. Observing 
the form of Friedmann equations ^, ([S]), as well as of 
P^ . we deduce that such a behavior can be easily ob- 
tained in principle. In the following two subsections we 
proceed to the detailed investigation in the case of scalar 
and /(i?)-modified cosmology. 



A. Scalar cosmology with Lagrange multiplier 

In order to provide a simple realization of cyclicity in 
this scenario, we start by imposing a desirable form of 
H{t) that corresponds to a cyclic behavior. We consider 
a specific, simple, but quite general example, namely we 
assume a cyclic universe with an oscillatory scale factor 
of the form 



a{t) = A sin(ajt) 4- Oc 



(19) 



where we have shifted t in order to eliminate a possible 
additional parameter standing for the phase^. Further- 
more, the non-zero constant Oc is inserted in order to 
eliminate any possible singularities from the model. In 
such a scenario t varies between — oo and ~foo, and t = 
is just a specific moment without any particular phys- 
ical meaning. Finally, note that the bounce occurs at 
<iB{t) = ttc — A. Straightforwardly we find: 






Au! cos{ujt) 
Asiii{ujt) + Uc 
Alu'^ [A + flc sin(a;i)] 
[Asin(ci;<) + ad 



(20) 
(21) 



Concerning the matter content of the universe we assume 
it to be dust, namely Wm ~ 0, which inserted in the 
evolution-equation pm + 3-ff(l -f- Wjn)pm = gives the 
usual dust-evolution p„i = pmo/a'^, with p„iQ its value at 
present. 

Now, we first consider (^ to be a canonical field, that 
is we choose aj((^) = 1. Concerning the potential C/(0) 
that will give as the solution for 0(t) according to ([H]), we 
choose a simple and easily-handled form, namely 



[/(0) 



m 



(22) 



where m is a constant with mass-dimension. In this case, 
(13) leads to 



also the A(i) evolution, which according to (|lip reads 



A(i) 



-1- 



Au'^ [A + Gc sin (wi)] 
[j4sin (cut) + flc] 



1 _3 

— jPmo{Asin{ujt) + ac} . 



(25) 



In order to provide a more transparent picture of the 
obtained cosmological behavior, in Fig. [T] we present the 
evolution of the oscillatory scale factor (fT9|) and of the 
Hubble parameter (|20|) . with A = 1, ui = 1 and Uc = 1-3, 
where all quantities are measured in units with SttG — 1. 
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rri^t, 



(23) 



having chosen the -I- sign in ©. As we have said in 
subsection III Al such a simple relation for 0(i) allows to 
replace t by (f)/m'^ in all the aforementioned relations. 
Therefore, substitution of dlQl) , dH]) and ([22]) into (|T2| . 
and using instead of t, provides the corresponding ex- 
pression for V{(j)) that generates such an iJ(i)-solution: 



F(0) 



Alu^ 



A + acSinUo:;^] 



Ai 



Aujt 



- {-^) 



Asmiuj:^] + flc J 



(24) 



Note that in the case of dust matter, the reconstructed 
potential does not depend on the matter energy density 
and its evolution. Finally, for completeness, we present 



^ Note that if the average of an oscillatory scale factor keeps in- 
creasing throughout the evolution, it could yield a recurrent uni- 
verse which unifies the early time inflation and late time accel- 
eration (23. [2I. 



Figure 1: The evolution of the scale factor a{t) and of the 
Hubble parameter H{t) of the ansatz ifjffi) . with A= 1, u = 1 
and Qc ~ 1.3. All quantities are measured in units where 
87rG= 1. 
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Figure 2: V{(t>) from \24^ that reproduces the cosmological 
evolution of Fig. [^ withuj{(f)) = 1, U{cj>) = m*/2 andm = 1.2. 
All quantities are measured in units where SttG — 1. 



Note that H{t) by construction satisfies the requirements 
for cycUcity, described in the beginning of this section. 
Furthermore, in Fig. [2] we depict the corresponding po- 
tential ¥{(/)) given by ([M|) . From these figures we observe 
that an oscillating and singularity-free scale factor can be 
generated by an oscillatory form of the scalar potential 
V{(j)) (although of not a simple function as that of a{t), 
as can be seen by the slightly different form of V{(t>) in its 
minima and its maxima). This y(0)-form was more or 
less theoretically expected, since a non-oscillatory V{(j)) 
would be physically impossible to generate an infinitely 
oscillating scale factor and a universe with a form of time- 
symmetry. Note also that, having presented the basic 
mechanism, we can suitably choose the oscillation fre- 
quency and amplitude in order to get a realistic oscilla- 
tion period and scale factor for the universe. Finally, we 
stress that although we have presented the above spe- 
cific simple example, we can straightforwardly perform 
the described procedure imposing an arbitrary oscillat- 
ing ansatz for the scale factor. 

The aforementioned bottom to top approach was en- 
lightening about the form of the scalar potential that 
leads to a cyclic cosmological behavior. Therefore, one 
can perform the above procedure the other way around, 
starting from a specific oscillatory V{(l)) and resulting to 
an oscillatory a{t), following the steps described in the 
first part of subsection III Al As a specific example we 
consider the simple case 

V{(j))=Vosm{uJvc^) + Vc, (26) 

with [/((/>) chosen as in ([22)1 and thus ([25]) holds too. De- 



spite the simplicity of V{(j)), the differential equation PU)) 
cannot be solved analytically, but it can be easily handled 
numerically. In Fig. [3] we depict the corresponding solu- 
tion for H{t) (and thus for a{t)) under the ansatz ([25]). 
with Vo = 3, ioy = 1 and Vc ~ 3, with a(0) = 3.2 and 
H{0) = —0.7 (in units where SttG =1). As expected, we 
do obtain an oscillatory universe. 



B. /(i?)-gravity with Lagrange multiplier 

In order to provide a simple realization of cyclicity in 
this scenario, we start by imposing a desirable form of 
H{t) that corresponds to a cyclic behavior. Due to the 
complicated structure of the reconstruction procedure, 
we will choose an easier ansatz comparing to the previ- 
ous subsection, which allows for the extraction of analyt- 
ical results. In particular, having described the general 
requirements for an H{t) that gives rise to cyclic cosmol- 
ogy in the beginning of the present section, we choose 
Hit) (and not a{t)) to be straightaway a sinusoidal func- 
tion, that is 

H{t)^AHSm{LJHt), (27) 

which gives rise to a scale factor of the form 



a{t) = Aho cxp 



Ah coii{ujHt) 



LOH 



(28) 



CO 




Figure 3: The evolution of the scale factor a{t) and of the 
Hubble parameter H{t), for a scalar potential of the ansatz 
126\) . with Vb = 3, LOv = 1 and Vc = 3, a(0) = 3.2 and 
H{0) — —0.7. All quantities are measured in units where 
SttG = 1. 



which is oscillatory and always non-zero. Inserting (j27p 
and its derivative into (fT7| we obtain R{t) as 

R{t) = GAhojh cos{ojHt) + 12 A% sm^{ojHt), (29) 

a relation that can be easily inverted giving 



t{R) 



arccos 



UJH 



3u}H + ^3(48^^ - 4i? + 3w 



Hi 



12A 



H 



where we have kept the plus sign in the square root. Fi- 
nally, f2{R) can be reconstructed using (fT8|). leading to 



MR) = - -ujji [48 Al-iR- 



3cj 



H 



2R + 3ujjj + ujHJSi^SAj^ - 4i? + 3ujjf) 



(30) 



In summary, such an ansatz for /2(-R) produces the cyclic 
universe with scale factor ([25)) . Note that /2 {R) has a re- 
markably simple form, and this is an advantage of the 
scenario at hand, since in conventional /(i?)-gravity one 
needs very refined and complicated forms of f(R) in order 
to reconstruct a given cosmological evolution |26| . How- 
ever, as we stated in the beginning of subsection III Bl in 
the /(7?)-version of Lagrange- modified gravity one can- 
not incorporate the presence of matter in a convenient 
way that will allow for an analytical treatment. The ab- 
sence of matter evolution in a cyclic scenario is a disad- 
vantage, since we cannot reproduce the evolution epochs 
of the universe. Therefore it would be necessary to con- 
struct a formalism that would allow for such a matter in- 
clusion, similarly to the case of Lagrange-multiplier mod- 
ified scalar-field cosmology. Such a project is in prepara- 
tion. 



IV. STABILITY ANALYSIS 

A central issue in all cyclic models is the stability anal- 
ysis of the cosmological perturbations along with back- 
ground evolution. In Newtonian gauge the linear metric 
perturbation is given by 



ds^ = -(1 + 2$)di2 + a^{t){l - 2*)ffx2 



(31) 



where the variable <I> is the so-called Newtonian potential 
which describes the scalar part of metric perturbation. In 
the following we will study the evolution of the Newto- 
nian potential near the bouncing point, for the two above 
scenarios respectively. 



A. Scalar cosmology with Lagrange multiplier 



which is applicable in all scales. Thus, in this solution 
there are two modes D(x) and ^(a;), which are arbitrary 
functions of spatial coordinates, and their explicit forms 
can be determined by certain boundary conditions. 

Observing the second term of the right-hand-side of 
(j33p . one may concern that the Newtonian potential 
might be divergent on the bouncing point, when H = 0. 
Fortunately, this does not happen in a generic bounce 
model. As it was shown in phenomenological studies of 
generic bounce scenarios |27| , one can impose an approxi- 
mate parametrization of the Hubble parameter as a linear 
function of the cosmic time, that is H = ^t, nearby the 
bouncing point t^ = 0, with ^ being a positive constant. 
Doing so, in the neighborhood of a bounce in a specific 
cycle the Newtonian potential can be solved as 



An explicit analysis on the Newtonian potential of the 
scenario of scalar cosmology with Lagrange multiplier 
was performed in j21| . A very remarkable feature of this 
model is that the system possesses no wavelike modes and 
thus the sound speed of perturbations is identically zero 
in all backgrounds. By virtue of this particular property, 
it could be possible to control the dangerous growth of 
unstable modes of cosmological perturbations in the con- 
tracting phases of the cyclic universe. 

Assuming that there is no anisotropic stress in our 
model, then 'J = <&. Consequently, we only have one de- 
gree of freedom of cosmological perturbations. According 
to the analysis of [21|, one deduces that when the back- 
ground evolution is dominated by the scalar field </>, the 
perturbation equation of the Newtonian potential can be 
expressed as 



$f 



D 



l-\/>--(e^-l) 



b — e 2 
as 



(34) 

where as is the value of the scale factor at the bouncing 
point. 



Relation ([34| presents convergent behavior, and thus 
the perturbations are able to pass through the bounce 
smoothly and without any pathology. Therefore, the 
scenario at hand indeed provides a satisfactory way to 
realize a cyclic picture without instability on its pertur- 
bations. 
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$-0, 



(32) 



B. /(i?)-gravity with Lagrange multiplier 



where the prime denotes the derivative with respect to 
In a. Note that the dust-like sound speed Cs = has 
been applied to eliminate the gradient term. Equation 
([5^ yields the generic solution [21 1 



, , , H r da 



a 



(33) 



When we are dealing with /(i?)-gravity with Lagrange 
multiplier, we cannot use the simplifying relation ^P — 
$, even under the assumption of zero anisotropic stress. 
Therefore, in order to study the number of degrees of 
freedom in this case it is convenient to expand the action 
into quadratic order, which can be formally expressed as 



Si^) = 



Tl,RR 



SiR^ + h 



gti:^ 



5g'"'5iR + 52R 



2 



59^") XC2 



(35) 



where SiR and 52R denote the perturbed Ricci scalar at 
first and second order respectively. Moreover, C2 is the 
perturbed constraint equation at first order, which form 



is given by 



-Sg'^'d^.Rd^R + df.Rd^'SiR + J2,r5iR = 0. (36) 



In principle, one could worry since the above action 
involves two scalar modes and higher derivative terms, 
which could imply instabilities. However, this is not 
the case since the higher derivative terms can be fixed 
by the perturbed constraint equations. In addition, the 
vanishing of the {i,j) component of the perturbed Ein- 
stein equation allows to eliminate one degree of freedom. 
Therefore, there is still only one mode of metric pertur- 
bation which is able to propagate freely. 

In the following, we would like to focus on the prop- 
agation of the metric perturbation around the bouncing 
point, which is crucial to the stability analysis of a cyclic 
scenario. For calculation convenience we assume that 
the bounce takes place slowly, and then the universe ap- 
proaches a static phase around the bounce asymptoti- 
cally. Under this assumption we can obtain the kinetic 
terms of the perturbation equation, which take the fol- 
lowing approximate form: 



Q-—Q 



1,R 



/: 



-Q 



= 0, 



(37) 



in which Q = $ -I- Vf. Although this equation is far from 
a complete form, one can still extract a few important 
issues. Namely, the sound speed of the perturbation is 
unity under the assumption of slow bounce, as it can 
be read from the coefficient before the gradient term. 
Furthermore, we confirm that there exist only a single 
degree of freedom in the ,f(R) cosmology with Lagrange 
multiplier. 

We mention here that we investigated the stability only 
under the assumption of slow bounce. And the result in 
this case illustrates that it is possible for the cosmolog- 
ical perturbation to evolve through the bounces. How- 
ever, it is necessary to point out that a more complete 
and generic analysis should be performed in order to con- 
strain the parameter space of the scenario. Such a general 
analysis is left for future investigation. 



leads to a rich cosmological behavior. In particular, one 
can obtain an arbitrary cyclic evolution for the scale fac- 
tor, by reconstructing suitably the scalar potential. As 
expected, an oscillatory scale factor is induced by an os- 
cillatory potential, and we were able to perform such a 
procedure starting either from the scale factor or from 
the potential. Additionally, the matter sector can be also 
incorporated easily, and this is an advantage of the sce- 
nario, since it allows for the successful reproduction of 
the thermal history of the universe. 

In the case of /(i?)-gravity, we considered a Lagrange 
multiplier field and a second form f2{R)- This sce- 
nario leads also to a rich behavior, and one can acquire 
cyclic cosmology by suitably reconstructing f2{R), which 
proves to be of a very simple form, contrary to the con- 
ventional /(i?)-gravity. However, the complexity of the 
model does not allow for the easy incorporation of the 
matter sector, since one cannot extract analytical results, 
and this is a disadvantage of the construction. 

In summary, we saw that Lagrange-multiplier modi- 
fied gravity may lead to cyclic behavior very easily, and 
the scenario is much more realistic in the case of scalar 
cosmology. Since a necessary test of every cosmologi- 
cal scenario is to examine the evolution of perturbations 
through the bounce [23|, we extended our analysis be- 
yond the background level, in both scenarios we consid- 
ered. For the case of scalar cosmology the perturbation 
behaves as a frozen mode without oscillations, since its 
sound speed is vanishing. For the case of /(i?)-gravity 
and under the assumption of slow bounce, we obtain a 
dynamical degree of freedom of which the sound speed is 
almost unity. Thus, in conclusion, it is possible for the 
cosmological perturbation to evolve through the bounces. 
However, a more generic perturbation analysis, beyond 
the slow bounce assumption, is needed, but since it lies 
beyond the scope of this work it is left for future investi- 
gation. 



CONCLUSIONS 



In this work we investigated cyclic evolutions in a uni- 
verse governed by Lagrange-multiplier modified gravity. 
In order to be more general we considered the scenario 
of modified gravity through Lagrange multipliers in both 
scalar-field cosmology, as well as in f{R) one. 

In the case of scalar cosmology, the use of a second 
field that is the Lagrange multiplier of the usual scalar. 
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